Magnetic order and moment distribution in doped spin-chain systems 
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We consider the effect of doping on the magnetic order in quasi one- dimensional (ID) antiferromagnets. First 
the detailed magnetic response of finite spin-1/2 chains in the presence of a staggered field is determined, which 
by itself shows an interesting crossover as function of length and field strength. We can then understand the 
ordering by including an effective coupling between the chains, which results in a Neel ordered phase at low 
temperatures, but with a broad distribution of magnetic moments that has to be determined self-consistently. 
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Doping and impurity effects in low-dimensional 
antiferromagnets remain of strong interest in the 
condensed matter physics community, spurred by 
high-T c superconductivity and other exotic ef- 
fects in those systems. Typically, antiferromag- 
netic three dimensional Neel order dos exist at 
extremely low temperatures, but is easily de- 
stroyed in such systems by a combination of quan- 
tum fluctuations, temperature fluctuations, and 
disorder effects. In this paper, we study how 
static impurities affect the Neel order in quasi- 
1D antiferromagnets, which are modelled by spin- 
1/2 chains with a weak coupling J' -C J to 
the neighboring chains 8 (before doping) H = 

In uncoupled chains it is known that long-range 
antiferromagnetic order is destroyed by quantum 
fluctuations even at T = 0, which results in anti- 
ferromagnetic correlations that decay with a pow- 
crlaw. However, a staggered field B can easily 
induce alternating order m a \ t — £ (5^) 
throughout the sample since the staggered sus- 
ceptibility diverges at T = pQ. 

In the case of quasi-lD antiferromagnets the 
coupling to the neighboring chains plays the role 
of an effective staggered field from the antiferro- 
magnetic background B = zJ'm a i t , which always 
implies an ordered state at T = p^. Here, z is 



the number of neighboring chains. Note, that the 
effective field neglects quantum fluctuations be- 
tween neighboring chains, but within each chain 
we are able to determine the alternating moment 
by a full quantum mechanical solution m a \t{B). 
We can then study the ordered state by solving 
the effective field equation self-consistently for a 
non-zero solution m 

m = m a i t (zJ'm). (1) 

For the undoped case it is known that the al- 
ternating magnetization is given by a scaling law 
TTiait(B) = aB 1 / 3 up to logarithmic corrections 



which results in m 



/ zJ' |3] in good 
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agreement with experiments 4 . 

In a doped system the impurities effectively cut 
the chain at low temperatures [H]. It is therefore 
necessary to calculate the alternating moment for 
finite chains and average over the disorder in or- 
der to arrive at a self-consistent distribution of 
moments m in Eq. JQ). Previous results have 
shown that finite chains have a different stag- 
gered susceptibility as a function of T depending 
if they contain an even or odd number of spins [H], 
which can be traced to the contribution of the 
zero- modes in the correlation functions 0. More- 
over, it has been found that the response is lower 
for finite chains which yields a distinct drop of the 
the Neel temperature with increased doping 6 in 
agreement with existing experiments 0] . The un- 
doped value of Tm scales roughly with J 1 £Q up 
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Figure 1. The alternating magnetization as a 
function of staggered magnetic field B and length 
L from numerics for open chains with L < 27 and 
B < 0.5 J. Inset: The corresponding distribution 
of moments in units of hgns according to Eq. JTJ 
for J' = 0.0008J. The non-uniform distribution 
of m along the chains was taken into account by 
including the known behavior near the edges |6I8| . 



to logarithmic corrections which we have quanti- 
tatively determined to give 0.6 < T^/zJ' < 1 in 
the range 0.1 > J'/J > 0.00001. 

The alternating moment in a finite field at 
T = in Fig. 1 shows again a distinct difference 
between even and odd chains and obeys a scaling 
behavior of the form 

malt (B) = L-^f(L^B^), (2) 

in agreement with a dimensional scaling analysis. 
Assuming a sharp distribution of lengths L ss l/p 
Eq. (J2J) also implies a scaling form for the average 
alternating moment 

m = L- 1/2 g(zJ'L). (3) 

The low doping limit hin^^oo f(x) — > ax re- 
covers the B 1 / 3 scaling behavior with a length- 
independent solution for m = a 2 / 3 \/zJ', which 
we see as a sharp peak in the moment distribution 
for low doping levels in the inset of Fig. 2] For 



larger doping p we expect a broader distribution, 
but the difference between even and odd chains 
also becomes important. The self-consistent so- 
lution to Eq. shows two distinct peaks cor- 
responding to the even and odd chain segments 
as displayed in the inset of Fig. ^ Interestingly 
the peak for the odd chains shifts to higher mo- 
ments (increased order parameter) as the doping 
increases, even though the Neel temperature is 
known to drop. In fact we know from Fig. ^ that 
m a it oc L -1 / 2 diverges for short odd chains. How- 
ever, in that limit quantum fluctuations between 
the shorter chains also become more important 
so it is not clear if this exotic "increased order 
from disorder" effect is robust or if the order pa- 
rameter always decreases with increasing doping. 
Qualitatively, however, the peak will indeed split 
for larger doping p > z J'/IO and the shape of the 
distributions in the inset is adequate, except for 
a doping dependent rescaling of m. 

In summary, we have studied the doping de- 
pendence of the Neel ordered phase at T = in 
quasi-lD antiferromagnets and showed explicitly 
how the order parameter is modified in the pres- 
ence of static disorder. The distribution of the 
magnitude of alternating moments broadens with 
doping and we find a characteristic double peak 
structure for p > z J'/IO. 
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